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ABSTRACT 



This report presents an algorithm and HP-67/97 programs 
for position determination with Loran-C chains. Operational 
data cards are prepared in advance for Loran-C triplets. 
Position determination is performed using a single program card 
and an appropriate operational data card. 



A. Introduction 



The Loran system is a radio aid to navigation which 
utilizes the principle of hyperbolic fixing. The locus of 
points for which the difference in arrival time of synchronized 
signals from a pair of transmitters is constant determines a 
hyperbolic line of positions (LOP) . The intersection of two 
hyperbolic lines of position from two pairs of transmitters 
determines position or a hyperbolic fix. That two pairs of 
stations are required for a fix does not necessarily mean that 
there are four separate stations, for one station of one pair 
may be colocated with one station of the other pair forming 
a Loran triplet (Figure 1). Triplets may be joined "end-to-end" 
by station colocation to form a Loran chain (Figure 2) . Loran 
chains are common on both the East and West Coasts of the 
North American continent. 

The early "Standard Loran" or Loran-A" operating at a 
frequency just below 2MHz is still in use in the Pacific area. 
The present day "Loran-C" operates at 100-kHz and is in use 
in both the Atlantic and Pacific Areas. The computational 
algorithm and programs described herein can be used for position 
determination with Loran-C triplets. Further information on 
the history, development and operation of the Loran systems 
may be found in References 1 and 2. 
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(a) Colocated Master Stations 



S 



1 



,S 



2 




(b) Colocated Slave Stations 




(c) Colocated Master and Slave 



Figure 1. Loran Triplets. 
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Figure 2. Loran Chain of Five Loran Triplets. 
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B. Program Description 



One program card and one operational data card 
(described below) are all that is required for on-location 
position determination from Loran triplet time-difference 
measurements. Two program cards are required to prepare 
operational data cards; these operational data cards should 
be prepared and validated prior to on-location navigational use. 
Thus although three program cards are described only one program 
card is required for navigation; two program cards are used to 
prepare operational data cards during or prior to mission planninc 
The function of each program card and its intended use follows. 



Program Card 1 . This program card is used to prepare 
master data cards. A master data card requires the following 
information for a master(M) station/slave (S) station pair: 

1. A M/S pair identification number. 

2. The quantity At which is the sum of the coding delay plus 
the one way base line time in microseconds. 

3. The latitude and longitude of the master station. 

4. The latitude and longitude of the slave station. 

Some preprocessing of these data is performed before the master 
data card is generated. The data generated require only one side 
of an HP-67/97 magnetic card for each M/S pair, thus a second 
M/S pair may be placed on side 2 of the card (thus conserving 
cards) if desired. It is envisaged that a master data card 
will be prepared in advance for each M/S pair that might be 
received within an area of operation. 
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Program Card 2 . This program card is used to prepare 



an operational data card for every Loran triplet within an 
operational area. Each operational data card contains data 
merged from the master data cards which contain M/S pair 
information for each pair of the triplet. These merged data 
are validity checked, colocation of master or slave determined 
and encoded. 

The only inputs required for this program are the two 
master data cards that comprise the Loran triplet. It is 
possible to prepare and store operational data cards rather 
than master data cards. This may be desirable if there is no 
scarcity of cards and storage space, however the number of 
possible Loran triplets is considerably larger than the 
number of M/S pairs. 

Program Card 3 . This program card is used in conjunc- 
tion with an operational data card for position determination. 
Required input is the indicated time difference T for each 
M/S pair of the triplet. Output is the computed latitude and 
longitude of the fix. Note: Every Loran fix has two possible 

solutions. The unwanted solution can almost always be rejected 
by inspection, however, if the stations of the Loran triplet 
are nearly aligned then either solution may be valid even 
though only one solution should be consistent with the flight 
plan . 
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C. HP-67/97 Calculator Programs 



1. User Instructions 



CARD 1 







Input 




Output 


Step 


Instructions 


Data/Units 


Keys 


Data/Units 


1. 


Read in program card (both 
sides) 








2. 


Input a unique ID number for 
the Loran pair* 


ID 


f a 


ID 


3. 


Input the coding delay At 


At 


f c 


At 


4a. 


Input the master station latitude 




t 


— 


b. 


and longitude (CHS for West) ** 


X M 


A 


— 


5a. 


Input the slave station latitude 


$s 


t 


— 


b. 


and longitude (CHS for West) 


x s 


C 


— 


6a. 


Run 


None 


E 


crd 


b. 


Pass a blank data card through 
the card reader. 









Note: Loran pairs are coded on the navigation maps using 

designators such as 9930X, 9930Y, 99302 and 9930W. It is 
suggested that the ID's for these pairs be coded as 
9930.1, 9930.2, 9930.3 and 9930.4, respectively. However, 
any consistent scheme is acceptable. 

* 

The format for position data input is of the form: 

+DDD.MMSSFF, where 

DDD denotes degrees 
MMM denotes minutes 
SS denotes seconds 
FF dnotes hundredths of a second. 

The minus sign (-) denotes Southern latitudes or Western 
longitudes . 
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CARD 2 



r - - 

Step 


Instructions 


Input 

Data/Units 


Keys 


Output 

Data/Units 


1. 


Read in program card 2 








2a . 


Start 


— 


E 


9.00 


b. 


9.00 will flash in the display. 
Insert a master data card 
containing the first pair of 
a Loran triplet. 








c . 


9.00 will flash in the display 
once more. Insert a master 
data card containing the second 
pair of the Loran triplet. 






9.00 


d. 


If the data form a proper 
triplet, "crd" will appear in 
the display. 






crd 


e . 


Pass both sides of a blank card 
thru the card reader to produce 
the operational data card for 
the Loran triplet. 










Should "error" appear in the display, then the two master data 
cards do not compare to form a Loran triplet. Both the latitude 
and longitude of the colocated stations must be identical on 
both master data cards in order to successfully produce an 
operational data card. 


Label the A key position with the identification number of the 
first Loran pair (the pair inserted in Step 2b) and label the B 
key position with the identification number of the second Loran 
pair (from Step 2c). 
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CARD 3 



Step 


Instructions 


Input 

Data/Units 


Keys 


Output 

Data/Units 


1. 


Read in both sides of the program 
card 3 . 








2. 


Read both sides of the operational 
data card for the Loran triplet 
that you are receiving. 








3a. 


Set to compute Solution A. 




f a 


— 


b. 


Set to compute Solution B. 




f b 


— 


4. 


Input the observed time delay 
from the first Loran pair. 


T 


A 


— 


5. 


Input the observed time delay 
from the second Loran pair. 


T 


B 


— 


6. 


Compute fix Latitude 

Longitude 




E 

R/S 


Latitude 

Longitude 


7. j 


Repeat from Step 2 with a new 
operational data card or from 
Steps 3 or 4 as required. 
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2 . Sample Problem 



CARD 1 







Input 




Output 


Step 


Instructions 


Data/Units 


Keys 


Data/Units 




In this series of examples we 
will prepare and use data 
cards for the Loran-C pairs 
9930X and 9930Y. 








1 . 


Read in program card 1 
(both sides) 








2. 


Input the ID for 9930X 


9930.1 


f a 


9930.10 


3. 


Input the coding delay At 
for 9930X 


36389.66 


f c 


36389.66 


4a. 


Input the master station 
latitude 


34.034604 


t 




b. 


and longitude (CHS for West) 


-77.544676 


A 


— 


5a. 


Input the slave station 
latitude 


46.463218 


t 




b. 


and longitude (CHS for West) 


-53.102816 


C 


— 


6a. 


Compute 


None 


E 


crd 


b. 


Pass a blank data card through 
the card reader. Label the 
card 9930X MASTER 






— 


7. 


Input the ID for 9930 Y 


9930.2 


f a 


9930.2 


8. 


Input the coding delay At 
for 9930Y 


52541.31 


f c 


52541.31 


9a. 


Input the master station 
latitude 


34.034604 


1 




b. 


and longitude (CHS for West) 


-77.544676 


A 


— 


10a. 


Input the slave station 
latitude 


41.151193 


1 




b. 


and longitude (CHS for West) 


-69.583909 


C 


— 


11a. 


Compute 


None 


E 


crd 


b. 


Pass a blank data card (or 
the second side of the card 
used in Step 6b) through the 
card reader. Label the side 
9930Y MASTER 








12. 


These two cards will be used 
in the next example. 
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CARD 2 



Step 


Instructions 


Input 

Data/Units 


Keys 


Output 

Data/Units 


1. 


Read in program card 2. 








2a. 


Start 


— 


E 


9.00 


b. 


While 9.00 is flashing in the 
display, insert the MASTER 
data card for station 9930X 
into the card reader. 








c. 


When 9.00 starts flashing in 
the display again, insert 
the MASTER data card for 
station 9930Y into the card 
reader. 






9.00 


d. 


"Crd* will appear in the 
display. 


— 




crd 


e. 

j 


Pass both sides of a blank 
card through the card 
reader. 

Label this card 9930X/9930Y 
OPERATIONAL DATA CARD. Then 
label the A key position 9990X 
and the B key position 9930Y. 
This card will be used in the 
next example. 






0.00 
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CARD 3 



Step 



Instructions 



Input 

Data/Units 



Keys 



Output 

Data/Units 



You are receiving 9930X and 
9930Y and wish to obtain a 
fix. 



1. Read in program card 3 
(both sides) 

2. Read in the operational data 
card for the triplet 
9930X/9930Y (both sides) 



3. 



Select Solution A. 



4. The indicated time delay is 
49400 ys from 9930Y. 

Input the indicated time delay. 

5. The indicated time delay is 
28800 ys from 9930X. 

Input the indicated time delay. 



49400 



28800 



f a 



B 



A 



0.00 



0.00 



6 . 



Solution A: 42°44 , 57 n N Latitude 
4 1 ° 07 ' 32"W Longitude 



E 42.4457 

R/S -41.0732 



Solution B: 27°00* 07"S Latitude 
' 102 o 27'12 M E Longitude 

< Since you are navigating over 
the North Atlantic, Solution A 
is the desired fix. 



7. 



Repeat from Step 2 with a new 
operational data card or from 
Steps 3 or 4 as required. 
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Program Storage Allocations and Program Listings 





Card 1 . 










Registers : 










RO: 


ID 


SO: 


L 


RA: 




Rls 


At 


SI: 


T 


RB: 




R2 : 


2c 


S2 : 


U 


RC: 


0 

m 


R3: 




S3: 


V 


RD: 


A6 

] 


R4 : 


0 M 

M 


S4: 


X 


RE: 


AX 


R5: 


X M 

M 


S5 : 


y 


RI : 


AX 


R6 : 


"MS 


S6: 


{ i d 






R7 : 


6 s 


S7 : 


ax' 

m 






R8 : 


A s 


S8 : 


d 






R9 : 


C SM 


S9 : 


f 







Initial Flag Status and Use: 



0: OFF, Unused 


2: 


OFF, Unused 


1: OFF, Unused 


3: 


OFF, Unused 


Display Status: 






DSP 4, FIX, DEG. 






User Control Keys: 






A: i t X 

M M 


a : 


Station ID 


B: 


b : 




C: d» s + s 


c ; 


At 


D: 


d: 




E: Prepare data card 


e : 
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Card 2 . 



Registers : 




R0 : 


+ ID 1 


SO : 


Rl : 


At^ 


SI: 


R2: 


2c-^ 


S2 : 


R3 : 




S3: 


R4 : 


CD 

> 

H 


S4 : 


R5: 


a ai 


S5 : 


R6 : 


^A1 


S6 : 


R7 : 


°B1 


SI: 


R8 : 


A B1 


S8 : 


R9 : 


^Bl 


S9 : 



Initial Flag Status and Use : 
0: OFF, Vertex determination 

1: OFF, Unused 



± id 2 


RA: 




At 2 


RB: 




2C 2 


RC: 






RD: 




6 A2 


RE: 


a = 21295.87 
P 


A A2 


RI: 


f = 1/298.26 


^A2 






0 B2 






A B2 






r B2 







2: 


OFF, 


Validity checking 


3: 


OFF, 


Unused 



Display Status : 

DSP 2, FIX, DEG 



User Control Keys: 



A: a : 

B: b: 

C : c : 

D: d: 

E: Run e : 
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Card 3 



Registers : 














RO : 


ID 1 




SO : 


ID 2 


RA: 


M 




Rl: 


At 1 




SI: 


At 2 


RB: 


u, 


N 


R2: 


(2c) 1 




S2: 


(2c) 2 


RC: 


D, 


d 


R3: 


A l' °2 ' 


P 


S3: 


A 2 


RD: 


A a 




R4 : 


CD 

h- 1 

II 

CD 




S 4 : 


9 2 


RE: 


a 

P 


= 21295.87 


R5 : 


A i - a f 




S 5 : 


A 2 


RF: 


f = 


1/298.26 


R6 : 


€ 1 




S6 : 


^2 








R7 : 


c l' c l' 


H 


S7 : 


C 2 








R8 : 


B 1 , S/a 


= r 


S8 : 


B 2 








R9 : 


a i 




S9: 











Initial Flag Status and Use: 



0: 


OFF, 


Soln A, Soln B 


2: 


OFF, 


M/S Vertex Flag 


1: 


OFF, 


Unused 


3: 


OFF, 


Unused 



Display Status : 

DSP 2, FIX, DEG 



a: Soln A 

b: Soln B 

c : 
d: 



User Control Keys : 
A: 

B 
C 
D 



E : Run 



e : 
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15 



Card 1; Prepare Master Data Card 
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Card 1: Continued. 
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Card 1: Continued. 
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Card 2: Prepare Operational Data Card 
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Card 2: Continued. 
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Card 3: Loran Fixing Program 
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Card 3 : Continued 
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Card 3: Continued 



D . Loran-C Fixing Algorithms 

The development of the Loran fixing algorithms in this 
report is presented in more detail in a companion report 
[Ref. 3] and will not be repeated here. 

The basic Loran-C equation [Ref. 4] can be written as 



T = 


[T s + P(T S )] - [T m + p(T M )] + [T b + p(T N )] + 6 (1) 



where 



T 


is the "indicated time difference" in microseconds. 


t m ,t s 


is the distance, in microseconds, from the master 
and the slave to the receiver, respectively, 


t b 


is the distance, in microseconds, between the 
master and the slave. 


6 


is the assigned coding delay, in microseconds, and 


p (T) 


is the secondary phase correction, in microseconds, 
for an all sea water path of length T. 



The quantity 

A t = [T b + p(T fi ) ] + 6 

is a constant for each master/slave pair. The following 
World Geodetic System 1972 (WGS 72) values have been adopted 
for Loran-C navigation [Ref. 4]: 
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Vq = 299792458 meters/second is the velocity of light 
in free space, 

n = 1.000338 is the index of refraction of the surface 
of the earth for standard atmosphere and 100kHz 
electromagentic waves, 

a g = 6378135.00 meters is the equatorial radius of the 
earth, and 

f = 1/298.26 is the flattening factor (1 - b/a^, where 
b is the polar radius) of the earth. 

Accurate formulas for computing the secondary phase 
correction p(T) are contained in Reference 4, but for use 
with the handheld calculator the following linear approximation 
[Ref. 3] will be used: 



where 

and 



P (T) = a 1 + a 2 T , 



a ± = -0.321, 
a 2 = 0.000635. 



Using this approximation, it is possible to solve Equation 1 
for the quantity Tg - T^. We find 



T s - t m = (T - At) / ( 1 + a 2 ) 



( 2 ) 



On the surface of a sphere a hyperbolic line of position 
can be represented by the equation [Ref. 3, page 175] 
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tan r 



cos 2a - cos 2c 

sin 2c cos ai + <; sin 2a 



(3) 



where the origin of the coordinate system is at the prime focus 
of the spherical hyperbola, 2c is the spherical arc joining 
the foci, 2a is a constant for any one LOP, and r and w 
are the spherical coordinates of a point on the LOP. If the 
base line of the coordinate system is the arc joining the foci 
then co is the spherical polar angle from the base line to a 
point P on the LOP and r is the spherical polar distance 
(or arc) from the prime focus to P. Using the Loran system 
we take c = +1 if the prime focus is at a master station 
and we take £ = -1 if the prime focus is at a slave station. 

If we take v = v^/h to be the velocity of 100kHz electro- 
magnetic radiation of the earth's surface then 



2a = v (Tg 



V //a e 



i 



or , using Eq . ( 2 ) , 

2a = (T - At) /a , (4) 

P 

where 

a (1 + a ) 

a = — 7 — = 21295.87 ys . 

P v 0 /n 

The baseline between master and slave can be obtained from 



2c = v T D /a . (5) 

B e 

Here 2c is computed by program card 1 (preparation of master 
data cards) using the algorithm in Section E. 
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Consider a Loran-C triplet with master stations colocated. 
Let and ? 2 denote the azimuth angles of slave 1 (S^) 

and slave 2 (S 2 )/ respectively, measured from North toward the 
East from the master stations (M) (see Fig. 3). Further, let 
a and r denote the azimuth and spherical polar arc (distance) 
of the receiver (R) from M. For this geometry, Eq. (3) can 
be written as 



B . 





tan r . = 
1 


-L 

C . cos (a-£ . ) + A . 

1 1 X 




(6) 


where 


A . = C . sin 2a . 
11 1 










B. = cos 2a. - 
1 X 


cos 2c. 

1 






and 


C . = sin 2c. 
1 1 








for the 


. th 

1 — Loran pair. 


i = 1,2. Since r = 


r l r 2 ' 


tan r . 
1 


can be 


eliminated in Eq 


. (6). The resulting 


equation 


can be 


rewritten as 










C cos a + 


S sin a = K , 




(7) 



where 





c = Bl c 2 


cos ^2 


- B 2 C 1 


COS 


i' 




s = Bl c 2 


sin C 2 


- B 2 C 1 


sin 


r i' 


and 


K = B 2 A 1 


- b x a 2 


• 
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N 




Figure 3. Geometry of a Loran Triplet and a Receiver. 
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If we define p > 0 and y by the equations 



p cos y = C, 

and (8) 

p sin y = S , 

then 

P = /C 2 + S 2 , 

and 

Y = qatn(S,C) . 



Here the function qatn(y,x) is the arctangent of y/x 
adjusted for the proper quadrant according to the signs of x 
and y. A compact form of this function is 



qatn(y,x) 



tan 



-1 



X 



X + 



-9 

10 t 



+ TTt (X < 0?) 



(x = 0?) 



where 



and 



t ( z) 



t ( z) 



1 when z is true 



0 when z is false. 



When convenient we will use the notation qatn(y/x) inter- 
changeably with qatn(y,x). The qatn function is equivalent 
to the polar angle obtained using the rectangular to polar 
conversion function on the HP-67/97. 

Now substitute Eq . (8) into Eq . (7) and solve for 

a = y + cos (k/p) (9) 
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to obtain the azimuth angle a of the two points of intersection 
of the LOP's. Finally we obtain a value for r by substituting 
each a into Eq. (5) . We find that 



r 



qatn 



B . 

1 

C . cos (a-£ . ) + A . 
1 l i 



for i = 1 or 2 . 



The distance and azimuth from M or the triplet vertex can be 
converted into the latitude and longitude of the two possible 
positions of R. 

The fixing algorithm then uses a and r in the direct 
solution algorithm of spheroidal geodesy (Section F) . 
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E. 



The Reverse (Inverse) Solution Algorithm 



This reverse solution algorithm is a modification of the 
first order in flattening (f) algorithm given by Thomas [Ref. 5, 
pp. 8-10]. Thomas' notation has been followed as closely as 
possible for ease of comparison of the algorithms. The qatn 
function is defined in Section D. West longitudes (X) and 
South latitudes ( 0 ) are negative. We are given the points 
P 2 ° n s Pheroid and are to find the 

distance S between the points and the forward and back 
azimuths, c* 12 and a 2 i* Given quantities are 4> 1 , A^ <f> 2 
and 2 . No assumptions about the relative location of P^ 
and P 2 are required. The modified reverse solution algorithm 
is : 



0^ = tan 1 ((l-f) tan <Jk] , i = 1, 2, 



0 = 

m 



(e +e_)/2, Ae = (e_-e,)/2, ax = a„ - a,, 

12 m 2 1 2 1 



2 2 2 2 

AA = AA/2, H = cos A0 - sin 0 = cos 0 - sin A0 = cos 0 cos 0 , 

m m m m m 1 2 

2 2 2 2 
L = sin A0 + H sin AA = sin (d/2) , 1 - L = cos (d/2) , 
m m 



-1 2 2 
d = cos (1 - 2L) , U = 2 sin 6 cos A0 /(I - L) , 

m m 

2 2 

V = 2 sin A0 cos 0 /L, X=U+V, Y = U - V, 
m m 



T = d/sin d , 6-d = f(TX - Y)/4, S = a (T-6-,d) sin d, 

i e 1 

F = 2 [Y - (1 - 2L) (4 - X) ] , G = fT/2 , 

Q = ~(FG tan AA)/4, AA 1 = (AA + Q) /2 
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t_ = qatn(-sin A0 cos AA ' , cos 0 sin AA 1 ) , 
1 m m m m 



t = qatn(cos A0 cos AA 1 , sin 0 sin AA 1 ) # 

2 m m m m 



a i2 = fc l + V “21 = fc l - fc 2 * 



This reverse solution algorithm is used by program card 1 

(preparation of master data cards) to compute the baseline 

distance 2c and the azimuths and between the 

MS SM 

master and slave stations of a Loran pair. 

Details of the modifications made to Thomas' algorithm 
are contained in Reference 3. 
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F. 



The Direct Solution Algorithm 



This direct solution algorithm is a modification of 
the first order in flattening (f) algorithm given by Thomas 
[Ref. 5, pp. 7-8]. Thomas' notation has been followed as 
closely as possible for ease of comparison of the algorithms. 
The qatn function is defined in Section D. West longitudes 
and South latitudes are negative. We are given the point 
Pl(4> ,A ) on the spherioid, where A are the geodetic 

latitude and longtiude (geographic coordinates); the forward 
azimuth a^ 2 and the distance S to a second point ^2^2'^2^ 
and from these we are to find the geographic coordinates 

\ 2 and the back azimuth a 2 i ’ The given quantities are 

, A a 1 _ and S. No assumptions about the relative loca- 
1 1 12 

tion of and P 2 are required. The modified direct 

solution algorithm is: 

5^ = tan ^ I ( 1 — f ) tan cjs^], M = cos 0^ sin 

2 

N = cos ^ cos a , c i = c 2 = ~ M ) /4 / 

D = 1 - 2c ^ - c^M, P = c^/D, = qatn(N, sin 0^) 

d = S/(a D) , u = 2 (G -d) , W = 1 - 2P cos u, 
e 1 

V = cos ( u + d) , Y = 2PVW sin d, Acs = d-Y , 

1 ^ j = qatn[-M, - (N cos A a - sin 0 sin Acs ) ] , 

2 1/2 

K = ( 1 — f ) |M + (N cos Ao - sin 0 sin Acs ) ] , 

‘P ^ = tan " | (sin 0 ^ cos Acs + N sin Acs ) /K] , 

• n = qatn (sin Acs sin a . cos 0, cos Acs - s in 0 sin Ao C os a ) 

12 1 1 12 ' 

H = CjAo. AA = An - H, A 2 =\ 1 + AA. 
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This direct solution algorithm is used by program 
card 3 (improved fix program) to compute the latitude and 
longitude of the receiver using the azimuth and range of the 
receiver from the Loran triplet vertex. 

Details of the modifications made to Thomas' algorithm 
are contained in Reference 3. 
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G. 



Discussion and Some Typical Results 



The HP-67 program design specifications of 
COMPATWINGSPAC [Ref. 6] are contained in the following 
statement . 

"There is a need for an HP-67 program that will compute a 
geographical position from two Loran delay rate readings. 

Several methodologies are available to compute the 
desired position but computational complexities increase 
with the desired accuracy and flexibility. The most 
desirable accuracy would be an error of less than 4 n.mi. 
at a range of 500 n.mi. with less error closer to the 
stations. It is likely that program length considerations 
will require that the station pairs have a common site 
(i.e. two slaves or two masters at the same location). 

This is not an unusual situation as evidenced by strings 
of station pairs along coast lines. A data card will 
probably be necessary for the station pairs to be used. 

However, more than one program card is unacceptable due 
to the decrease in functional utility when compared to 
the manual plotting method. As a final requirement, the 
fix should be obtainable on either side of the baselines 
onnecting the stations, and not limited to a geometric 
position relative to one side or the other of the 
stations . " 

It was further stated that the maximum computation time to 
obtain a fix be 1.5 minutes. 

It is felt that these design goals have been satisfied. 
Although one program is required to prepare master data cards 
for all Loran-C pairs and a second card is required to prepare 
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operational data cards, one each for every triplet, this pre- 
paration should be done only once. The data cards should be 
supplied to users verified and labeled, by the Fleet Mission 
Program Library. One program card and an appropriate oper- 
ational data card are all that is required for the fixing 
algorithm. 

The fixing algorithm will display one of the two 
possible receiver positions in 38 seconds following the entry 
of the time delay readings. Since there are situations in which 
either of the two solutions could be the valid solution; the 
decision of which solution to use should be left to the operator, 
not the program designer. 

Testing of the algorithm for all Loran-C triplets and 
positions relative to those triplets was too extensive a 
program to be carried out in the available time. Some "typical" 
scenarios however are presented in Tables I through IV. As can 
be seen all errors are all well within the design specifications 
of 4 n.mi at 500 n.mi range from the stations. The time delay 
values in these Tables were generated using a program discussed 
in Reference 3. It is recommended that the P-3 community test 
the algorithm for accuracy in known areas of operation and 
examine the results for possible regions in which the algorithm 
may fall outside the design requirements. Such testing should 
be compatible with the known "unreliable regions" shown on the 
Loran-C charts. 
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Table I. Moffett Field South 



Position 
Lat Long 


Indicated 

9940X 


Time Delay 
9940Y 


Lat (N) 


Fix 

Long (W) 


Error 

n.mi 


24°N 


122 °W 


27726. 19 


40912.76 


2 3 °59 ' 55 " 


122 °00 * 01” 


0.08 


26 


122 


27715.97 


40998.39 


25°59 ' 57” 


122 °00 ' 01” 


0.05 


28 


122 


27702.41 


41117.84 


27°59 ' 59” 


122 °00 * 00" 


0.02 


30 


122 


27683.53 


41291.85 


29°59 ' 59" 


122 °00 ' 00 n 


0.02 


32 


122 


27655.47 


41555.46 


32 °00 ' 00" 


122"00 1 00" 


0.00 


34 


122 


27609.63 


41959.57 


34 °00 ' 00" 


122 °00 ' 00 " 


0.00 


36 


122 


27523.56 


42544.11 


36°00 ' 00" 


121°59 ’ 59" 


0.01 


38 


122 


27334.61 


43248.22 


u> 

00 

o 

o 

o 

o 

o 


121° 59 1 58” 


0.03 



Table II. Moffett Field West 



i — 

Position 
Lat Long 


Indicated 

9940Y 


Time Delay 
9940W 


Lat (N) 


Fix 
Long (W) 


Error 

n.mi 


37°N 


1 2 2 ° W 


42892.86 


16257.23 


36°59 ' 59" 


122 °00 '01" 


0.02 


37 


125 


43056.68 


15765.13 


37°00 ' 00" 


125°00 ' 00" 


0.00 


37 


128 


43137.78 


15327.12 


37°00 ' 00" 


128 °00 ' 00” 


0.00 


37 


131 


43191.10 


14970.77 


37°00 ' 00" 


131°00 ' 00” 


0.00 


37 


134 


43232.38 


14683.74 


37°00 ' 00" 


134°00 1 00" 


0.00 


37 


137 


43267.42 


14449.40 


37°00 ' 00" 


1 37 °00 ’ 00” 


0.00 


37 


140 


43298.80 


14254.02 


37°00 ' 00" 


140°00 ’ 01” 


0.01 


37 


14 3 


43327.85 


14087.43 


37 °00 1 01" 


142 °59 ' 59” 


0.02 
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Table III. Brunswick Northeast 



Position 
Lat Long 


Indicated 

7930Z 


Time Delay 
9930X 


Lat (N) 


Fix 

Long (W) 


Error 

n.mi 


60°N 


30°W 


52437.86 


28451. 72 


60°00 ' 03 " 


29°59 1 32" 


0.24 


58 


35 


51960.93 


28391. 50 


o 

o 

o 

o 

o 

00 

LO 


34 °59 ' 46" 


0.11 


56 


40 


50992.37 


28359. 15 


55°59 ' 59" 


39°59 1 54 " 


0.06 


54 


45 


49292.46 


28370. 85 


53°59 ' 59" 


44°59' 57" 


0.03 


52 


50 


47165.60 


28490. 64 


52 °00 ' 00" 


49°59' 59" 


0.01 


50 


55 


45236.59 


29070. 48 


o 

o 

o 

o 

o 

O 

LO 


55 °00 ' 00" 


0.00 


48 


60 


44505.60 


30991. 94 


48 °00 ' 00" 


60°00 ' 00" 


0.00 


46 


65 


44475.70 


33697. 14 


46°00 ' 00" 


65°00 ' 00" 


0.00 


44 


70 


44588.91 


36567. 42 


43 °59 ' 59" 


69°59’ 59" 


0.02 



Table IV. Jacksonville Southeast 



Position 
Lat Long 


Indicated 

9930W 


Time Delay 
9930X 


Lat (N) 


Fix 

Long (W) 


Error 

n.mi 


9 °N 


47 °W 


13058.04 


36466. 46 


8°59 1 


'19" 


46°59 ' 22" 


0.92 


12 


52 


12984.71 


37288. 35 


11 °59 1 


1 34" 


51 °59 ' 37" 


0.57 


15 


57 


12898.73 


38267. 58 


14 °59 1 


'44" 


56°59 ' 47" 


0.34 


18 


62 


12793.91 


39431. 32 


17 °59 1 


'52" 


61 °59 ' 54 " 


0.16 


21 


67 


12656.52 


40794. 36 


20 °59 1 


'56" 


66°59 ' 57" 


0.08 


24 


72 


12451.30 


42330. 55 


23°59 1 59" 


71059.59.. 


0.02 


27 


77 


12097.12 


43876. 62 


27 °00 ' 


01" 


0 

0 

0 

0 

0 


0.02 


30 


82 


12973.95 


44768. 53 


30°00 ' 


'01" 


82 °00 ' 06" 


0.09 
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APPENDIX. Loran-C Station Parameters 



The following list contains the pertinent parameters 
for each Loran-C station pair. This list was compiled from 
data in Reference 4. Each column contains the following 
information : 

1. The Loran-C station pair designator 

2. At, the sum of the coding delay plus one way baseline 
time, including the secondary phase correction for an 
all seawater path, in microseconds. 

3. The master station latitude. 

4. The master station longitude. 

5. The slave station latitude. 

6. The slave station longitude. 

In this list, negative longitudes are West and positive 

longitudes are East. If desired, this convention may be 

reversed since the algorithms are independent of such external 

conventions; if this is done, care should be taken that the 

signs of all longitudes in the list are reversed. In columns 

3 through 6 the latitudes and longitudes appear to be in 

decimal form, but the actual format is DDD.MMSSFF (which is 

compatible with the HP-67/97 H.MS input mode) where 

DDD designates degrees, 

MM designates minutes, 

SS designates seconds, and 

FF designates hundredths of seconds. 



39 



4390X . 
4990V, 
593 OX . 
5930V, 
5990X . 
5990V, 
5990Z, 
7930N , 
79307: 
79-J UZ . 
" ? 9d0X , 
7930V , 
797SW , 
7 97 OX , 
7 970V , 
7'-'70Z, 
73001-1 , 
7980X , 
703 OV 
T’^SOZ 
"T 1 399 , • 
" 330 V , 
, P ■' _ ; 

■ 

J - ;■ n< 1 



-■ ~i u.L M l j 



-! J. 



" 7' “-f! 

_ i-i 

' J'f 0;, . 

“ Zf 7 h 1 ' 

9?t 7C 
"'970W 

j -t ‘p “h 1 

"'970V ■ 
°?79Z.. 

-i :30X , 






15972. 
34253. 
13131 . 
23755. 
1 3343 . 
23927. 
42266 . 
15863. 
27303. 
43212. 
1 3304 . 
29651 . 
30065 . 
15043. 
43944. 
63216. 
12309. 



45201 . 
61542. 

^ •,-n C - lr 
1 u ) J L 



cr j_- q q q ^ 

14 - 55 . 

31162. 

1 —'—'CT 
~t i 1 JI ■ 



l"'f 95. 







i -■« -• i 




3 9 7 s £ 




"t* ’*7' 

iL. •— 1 


16. 


444395 .• 


- 1 69 . 


3031 


20, 


• 20 


. 144916, 


-155 


. 5309 


7£i 


13 


1 6 . 


444395 , 


-169. 


3031 


20 , 


"1 i“i 
1 ^ i^i 


.234177, 


-173 


• 1 r’ 30 


20 


IJ IJ 


46 4 


432720, 


-067. 


5537 


71, 


.41 


. 151193, 


-869 


. 583909 


02 


46 « 


432720, 


-067. 


r r -? 

_i _i .j, | 


71 , 


46 


. 463218, 


-053 


. 1023 


16 


60 


51. 


575878, 


— | Z'O 


2202 


24, 


. 55 


. 262085 , 


-131 


. 151965 


06 


51. 


cr -7 cr “? 

•J f JO 1 O .■ 


-122. 


2202 


24, 


.47 


. 034799, 


-119 


. 4439 


cr q 


63 


51. 


575378, 


-122. 


2202 


24, 


, 50 


. 362972, 


-127 


.2129 


35 


02 


59. 


591727, 


-045. 


1027 


47, 


■ 64 


. 542653, 


-023 


.5521 


75 


l i‘ 


59. 


591727, 


-045. 


1027 


47, 


1 •-.* >_ 


. 175963, 


-007 


. 0426 


71 


20 


59. 


591727, 


-045. 


1027 


47, 


' 46 


.463213, 


-053 


. 1 028 


1 6 


45 


60 - 


1 9423 1 , 


-142. 


4331 


90, 


.57 


. 26202 1 , 


-152 


.2211 


U* Cm 


14 


6 J . 


1 9423 1 , 


-142. 


4331 


90 , 


. 55 


. 262085 , 


-131 


. 151965 


fc*4 


62 . 


175963, 


-007. 


0426 


71 , 


. 54 


. 432930, 


+ 008 


. 1736 




10 


62. 


175968, 


-007 . 


y4z6 


71, 


1 hn 


. 3306 1 5 , 


+014 


.2747 


80 


cr “i 


62. 


175968, 


-007. 


0426 


71 


• 64 


c a y- cr 
- .' 


-023 


.5521 


75 


39 


£2 . 


175963, 


-007. 


0425 


71 . 


.70 


cr a cr -i 

- Jl-Jib 1 .• 


-008 


. 4353 


69 


54 


30. 


593374 , 


-085. 


1 @09 


30, 


■ 30 


. 433302 , 


-090 


. 4943 


60 


•“« i"i 


30. 


593374 , 


-035. 


1009 


30 , 


' 2'b 


.315501, 


-097 


. 5000109 


O O 


30. 


593374, 


-085. 


1 009 


30. 


• cl 1 


. 015849, 


-080 


. £16 5 3 


52 


79 


30. 


593374 , 


-035 . 


1009 


30, 


.34 


. 034604 , 


-077 


. 5446 


76 


37 


i”i 


522051 . 


0-16. 


4305 


y6 j 


•-.cr 
1 -O' ■_.! 


. 312033, 


012 


. 3129 


9 b 


30 


Z‘ C* . 


522061 , 


015. 


4305 


96 . 


. 48 


. 532095, 


027 


.5201 


52 


£Q 


\ 0 ■ 


52206 1 , 


016. 


4385 


96. 


.42 


. 033649, 


003 


. 1215 


98 


•J A 


'74 

... e 


510754, 


-087. 


2912 


14. 


. 30 


. 593374 , 


-035 


. 1009 


38 


06 


-;Q 

■j — • . 


510754, 


-037. 


2912 


14, 


■ 42 


. 425060, 


-076 


. 4933 


86 


74 


3 9 . 


510754, 


“ y y r* ■ 


2912 


14, 


.48 


. 364984 , 


-094 


.3318 


47 


51 


34 ! 


h y 4 b y 4 .1 


— 0 r* i*’ ■ 


5446 


7 6 , 


• c_ r 


.015849, 


-030 


. 0653 


52 


66 


34 . 


9 34604 ■ 


0 r r ■ 


5446 


“ 1 

1 ' b . 


.46 


. 463213, 


-053 


. 1 028 


1 6 


”1 -j 


34 . 


tio46a4 


— 0 r r • 


5446 


7 6 j 


.41 


. 151193, 


-069 


. 583909 


r .2 


34 . 


034604 . 


-077. 


5446 


— "i 

t b . 


. 39 


.510754, 


-037 


.2912 


14 


0 


*4 ^ 


330662 • 


-113. 


4956 


0 7 


. 47 


. 034799, 


-119 


. 4439 


53 


50 


— 1 — * c 


330662 .• 


-113. 


4956 


•Z i r j 


. 33 


. 465699, 


-122 


. 2944 


53 


jfci 


'-i 


ootlbbii .• 


-113. 


495b 


“;"7 


. 35 


. 191313, 


-114 


.4817 


4>! 


20 


42 ! 


425060 , 


— 0 7 6 . 


4933 


36 j 


. 46 


. 432720, 


-067 


cr cr -7 


71 


-4 q 


42. 


425060, 


-076. 


-1 933 


36 j 


41 


. 151193, 


-069 


. 533989 


65 


42. 


425060 ' 


-076. 


4933 


3 6 j 


.34 


. 034604, 


-077 


. 5446 


7 6 


06 


42. 


425060 , 


-076. 


4933 


36 j 


. 39 


.510754, 


-037 


.2912 


14 


94 


24. 


48041 , 


141 . 


19290 , 


.24 


. 17077 , 


153 


.5851 


5 


12 


24. 


43041 , 


141 . 


1 9290 j 


.42 


. 443700, 


143 


. 430906 


1 7' 

i. - 


24. 


43041 , 


141 . 


19290 , 


26 


■ 362433, 


123 


. 0856 


21 


73 


24. 


43641 , 


141 . 


19290 , 


. 09 


. 324566 j 


133 


. 0955 


•~l •*■*! 

«iO 


cl 


cr — ’ 


090933 , 


-17G. 


1459 


31, 


52 


.494505, 


+ 173 


. 1052 


31 


09 


cr — 1 

J 1 . 


03Q93 8 , 


-170. 


1459 


31 , 


65 


. 144012, 


— 1 66 


.531447 


10 


cr ”7 

» . 


096933 , 


- 1 70 . 


1459 


8 1 , 


cr -? 
. 1 


. 26202 1 , 


-152 


.221 1 


2 ci 



40 



Coverage of Loran-C Systems 



Station 


Location 


4990 


Central Pacific 


5930 


East Coast, Canada 


5990 


West Coast, Canada 


7930 


North Atlantic 


7960 


Gulf of Alaska 


7970 


Norwegian Sea 


7980 


Southeast U.S.A. 


7990 


Mediterranean Sea 


8970 


Great Lakes 


9930 


East Coast, U.S.A. 


9940 


West Coast, U.S.A. 


9960 


Northeast U.S.A. 


9970 


Northwest Pacific 


9990 


North Pacific 
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COMNAVSURFI.ANT 1 

Norfolk, FA 23511 
Attn: Code N612 

Deputy Commander 1 

Operational Test and Evaluation Force Pacific 

NAS North Island 

San Diego, CA 92135 

Attn: Code 202 

VP-24 1 

NAS 

Jacksonville, FL 32212 
Attn: LTJG Kenneth D. Walker 
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Assistant Wing TAC D&E Officer 1 

ASW Operations Department 

NAS 

Cecil Field, FL 32215 
Attn: Code VS-24 

Naval Electromagnetic Spectrum Center 1 

Naval Communications Unit 
Washington, D.C. 20390 
Attn: CDR Claude I aVarre 

VS-Support Unit 1 

NAS 

Cecil Field, FL 32215 
Attn: CWO Wentworth 

PROCAL-Code 18 1 

Fleet Readiness Office 

NOSC 

San Diego, CA 92152 
Attn: Jim Grant 

Chief of Naval Materiel 1 

MAT 08D1 

Department of the Navy 
Washington, D.C. 20360 
Attn: CDR Philip Harvey 

LT Jan Smith 1 

COMCRUDESGRU TWO 

FPO New York, NY 09501 

Dr. Thomas Burnett 1 

Dept, of Management Science 

OR Study Group 

U.S. Naval Academy 

Annapolis, MD 21401 

U. S. Naval Oceanographic Office 1 

Washington, D.C. 20309 
Attn: Paul D. Thomas 

Navigation Department 1 

Defense Mapping Agency 

Hydrographic/Topographic Center 

Washington, D.C. 20315 

Attn: G. R. Deyoung 
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LT Michael D. Redshaw 1 

Helicopter Antisubmarine Squadron 7 
FPO New York 09501 

Mr. Raymond F. Fish 1 

Naval Underwater Systems Center 
Newport, RI 02840 

LT Kenneth W. Peters 1 

Commanding Officer 

U. S. Naval Air Facility 

APO New York 09406 

LT Peter W. Marzluff 1 

Surface Warfare Officer's School 
Newport, RI 02840 

LT Ray C. Pilcher, Jr. 1 

Surface Warfare Officer's School 
Newport, RI 02840 

LT William M. Yerkes 1 



ASWOC NAF 
S^'gonella Sicily 
FPO New York 09523 

LT A. J. Koorey 1 

Commander 

Cruiser-Destroyer Group Two 
FPO New York 09501 

LCDR E. G. Schwier 1 

USS Estocin ( FFG 15) 

FPO New York 09501 

LCDP R. R. Hillyer 1 

Code N311 

Commander, Carrier Group Four 
FPO New York 09501 

LCDR J. W. Pattison 1 

Code N321 

Commander, Carrier Group Eight 
FPO New York 09501 

LT Richard C. Myers 1 

ASWOC NAF LAJES 

Azores 

APO New York 09523 
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LT R. N. Christianson 
ASW 

Naval Postgraduate School 
Monterey, CA 93940 

LT L. R. Erazo 
ASWOC SIGONELLA 
FPO New York 09523 

LT R. E. Springman 
USS Lexington (ANT-16) 

Naval Air Station 
Pensacola, FL 32504 

LT Philip D. Ward 
Patrol Squadron 26 
Special Project Det 
Naval Air Station 
Brunswick, ME 04011 

LT Robert J. Knight 

USS John F. Kennedy (CV-67) 

FPO New York 09501 

LT Richard 0. White, Jr. 

USS Dwight D. Eisenhower (CVN-69) 
FPO New York 09501 

LT David E. Hebdon 

USS John F. Kennedy (CU-67) 

FPO New York 09501 

LCDR Michael A. Santoro 
Air Antisubron 41 
NAS North Island 
San Diego, CA 92135 

LCDR Robert M. Hanson 
Helasron Four 
NAS North Island 
San Diego, CA 92135 

LT Michael D. Thomas 
2476 Ridgecrest Ave . 

Orange Park, FL 32073 

LT Carl E. Garrett, Jr. 

Class 65, Dept. Head Course 
Surface Warfare School, Naval Base 
Newport, RI 02840 
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Naval Postgraduate School 
Monterey, CA 93940 

Attn: Prof. R. N. Forrest, Code 55Fo 1 

Prof. R. H. Shudde, Code 55Su 200 

R. J. Stampfel, Code 55 1 

LT M. D. Clary 1 

USS PAIUTE ( ATF-159 ) 

FPO New York 09501 
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